arXiv:1501.06755vl [physics.flu-dyn] 27 Jan 2015 


Concentrations of inertial particles in the turbulent wake of an immobile 
sphere 

Holger Homann and Jeremie Bee 

Laboratoire J.-L. Lagrange UMR 7293, Universite de Nice-Sophia Antipolis, CNRS, Observatoire de la Cote d'Azur, 

Bd. de VObservatoire, 06300 Nice, France 

(Dated: 28 January 2015) 

Direct numerical simulations are used to study the interaction of a stream of small heavy inertial particles with 
the laminar and turbulent wakes of an immobile sphere facing an incompressible uniform inflow. Particles 
that do not collide with the obstacle but move past it, are found to form preferential concentrations both 
in the sphere boundary layer and in its wake. In the laminar case, the upstream diverging flow pattern is 
responsible for particle clustering on a cylinder that extends far downstream the sphere. The interior of this 
surface contains no particles and can be seen as a shadow of the large obstacle. Such concentration profiles 
are also present in the case of turbulent wakes but show a finite extension. The sphere shadow is followed 
by a region around the axis of symmetry where the concentration is higher than the average. It originates 
from a resonant centrifugal expulsion of particles from shed vortices. The consequence of this concentration 
mechanism on monodisperse inter-particle collisions is also briefly discussed. They are enhanced by both 
the increased concentration and the presence of large velocity differences between particles in the wake. 

PACS numbers: 52.30.-q, 52.65.-y, 52.30.Cv 


I. INTRODUCTION 

Industrial and environmental problems often require modeling the hydrodynamic interactions between particles sus¬ 
pended in a turbulent flow. Sometimes it is necessary to accurately model the collision efficiencies between droplets. 
Instances are rain triggering by droplet coalescences in warm clouds 1 1 or planet formation by dust accretion in circum- 
stellar disksP In such problems one encounters different regimes: For very tiny particles associated to small values 
of the Reynolds number, the flow surrounding the particles is purely viscous and evolves according to the Stokes 
equation. When dealing with a large number of particles, this approach leads to model the collective hydrodynamic 
interactions in terms of Stokesian dynamics ! 3 In turbulent carrier flows, one then generally approximates the flow as 
the superposition of the small-scale effects of particles and the inertial-range Navier-Stokes dynamics. 4 However, for 
larger particles with non-negligible Reynolds numbers, the boundary layer and wake can become rather complicated, 
so that for example even the drag force experienced by such particles is up to now only heuristically modeled. There¬ 
fore very little is known on the challenging problem of modeling the hydrodynamic interactions between two particles 
with Reynolds numbers implying that at least one particle has a turbulent wake. 

In this work we focus on the asymptotic case of particles with very different sizes and investigate how a fixed 
sphere developing either a laminar or turbulent wake interacts with a stream of very small heavy particles. Inertia 
makes the trajectories of the small heavy particles deviate from fluid element trajectories and this possibly leads to 
particle-sphere collisions. One question is then to understand the rate at which such collisions occur and how it is 
accelerated or depleted as a function of the obstacle Reynolds number. This problem has important applications in the 
determination of dust accretion by planetesimals in early stellar systems or in the wet deposition of tiny aerosols by 
raindrops and is studied elsewhere. In the present work, we are concerned with the dynamics of those particles which 
pass the obstacle and interact with its wake. Using direct numerical simulations, we study their spatial distribution 
and dynamics downstream the spherical obstacle . 

It is well known, that inertial particles do not exactly follow the trajectories of fluid elements and distribute with 
an inhomogeneous density. 5 7 Heavy particles are ejected from coherent rotating structures of the flow by centrifugal 
forces and concentrate in stretching regions. The creation of concentrations and voids is correlated to the local struc¬ 
ture of the fluid velocity gradient tensor; particles concentrate in regions where it has real eigenvalues and voids appear 
in places where two eigenvalues are complexPFigure[l]shows a snapshot from one of our simulations displaying small 
inertial particles in the turbulent wake of a sphere. Detached coherent vortices that are advected downstream by the 
mean flow eject particles and create voids. At the border of these voids, clumps of particles can be observed. These 
clusters trigger collisions that might be important in applications. Cloud water droplets can be this way concentrated 
by a falling raindrop, eventually coalesce, grow in size and influence the rain formation process. 

Tang et al. 9 and Yang et al ) 10 studied numerically and experimentally the influence of advected large scale structures 
on particle dispersion in plane mixing layers and wakes. They show that the development of highly organized patterns 
of inertial particles is determined by the evolution and interaction of these structures. Via the measurement of the 


2 



FIG. 1. Snapshot of the position of particles (with St = rUdd = 3.2, where XJ C is the homogeneous inflow velocity and d is 
the diameter of the spherical obstacle) in a small slice passing through the axis of symmetry. The sphere Reynolds number is here 
Re s = U c dlv = 400. The background shows the values of the discriminant A = (det[cr]/2) 2 - (tr[cr 2 ]/6) 3 which is positive in 
rotating regions and negative in stretching regions. 


fractal correlation dimension they found that particles with a certain mass are the most influenced. Our work is also 
related to an experimental study of Jacober and Matteson. 11 In a wind tunnel they measured average concentrations of 
inertial particles in the turbulent wake of a sphere with Reynolds numbers ranging from 23000 to 110000. Despite the 
fact that these Reynolds numbers are out of reach by direct numerical simulations we find good qualitative agreement 
with their data. 

This article is organized as follows. Section [II] gives a short description of the numerical method applied. Sec- 
i III is dedicated to the mean concentrations of small particles in the turbulent wake, while Section [TV] presents an 
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estimation of inter-particle collisions. Section |V| encompasses concluding remarks. 


II. THE NUMERICAL METHOD 

In order to conduct direct numerical simulations of a three-dimensional hydrodynamic flow around a spherical 
object with no-slip boundary conditions at its surface, we use a combination of a standard pseudo-Fourier-spectral 
method with a penalty method. 12 More precisely, we solve the incompressible Navier-Stokes equations 
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d t u + u • Wu =- Wp + vV I u + /, V • u = 0, (1) 

Pf 

for the fluid velocity u , where pf is the fluid density, v its kinematic viscosity and / a force maintaining a uniform 
and constant inflow. This mean flow is imposed by keeping constant the zero mode of the stream-wise component of 
the velocity; the mean velocity is denoted by Uc- Equation Q is associated to the no-slip boundary condition on the 
surface of the large particle at rest 


u(x, t) = 0 for \x-X s \< d/2. (2) 

Here X s denotes the position of the particle center and d its diameter. Homogeneous inflow conditions are achieved 
by using again the penalty method to remove the velocity fluctuations originating from the sphere wake at the exit of 
the computational domain, before they are injected upstream by periodic boundary conditions. For the time integration 
of 0 we use a third order Runge-Kutta scheme. The grid resolution is chosen to resolve all small scales: those of the 
boundary layer of the sphere and all turbulent scales in the wake. Details of the method and benchmarks can be found 
in Homann et al . 12 

For the small inertial particles, we consider spherical particles with a radius a much smaller than both the viscous 
boundary layer and the smallest active scale of the sphere turbulent wake (the local Kolmogorov dissipative scale 
rj = (y 3 /6) 1/4 , where 6 is the local energy dissipation rate). In this limit, particles can be approximated by point 
particles. Further, we assume that these particles move sufficiently slowly with respect to the fluid and that their mass 
density p p is much higher than the fluid density pf. With these assumptions the dominant hydrodynamic force exerted 
by the fluid is a Stokes viscous drag, which is proportional to the velocity difference between the particle and the fluid 
flow 1314 


X = (3) 

where the dots stand for time derivatives. The quantity r = 2p v a 2 / (9pfv) is called the particle response time and is 
a measure of the particle inertia. It is the typical relaxation time of the particle velocity to that of the fluid. We also 
assume that the particles are sufficiently diluted to neglect any interaction among them and any back-reaction on the 
flow. Usually, particle inertia is measured in terms of the Stokes number St = t / Tf defined by non-dimensionalizing 
their response time by a characteristic time scale Tf of the carrier flow. The different time scales involved here play 
a role depending on which aspect of the problem we are interested in. In most cases, we use Tf = d/Uc , which 
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corresponds to expressing the particle response time in units of the time needed to be swept over a distance d by 
the large-scale fluid velocity Uc- Particles with small St are closely correlated to the flow and are swept around the 
obstacle. Large St- particles will preferentially collide with it. Only when concerned with the particle dynamics in the 
turbulent wake we will use St^ = r / Tf with Tf = = (vie) 112 (the local Kolmogorov time scale). 

We report results of three simulations corresponding to three different values of the obstacle Reynolds number 
Re s = Uc d/v = 100, 400, and 1000. With this choice we cover various wake types. For the smallest Re s the wake 
is steady and laminar. The intermediate value is slightly above the onset of chaotic shedding of vortices and thus 
slightly turbulent. For Re s = 1000 the wake is fully turbulent. In all cases, we consider streams of heavy particles 
with response times r between 0.04 and 40.96, corresponding to Stokes numbers in the range 0.05 < St < 63. The 
main parameters of the simulations are summarized in Tab. |T| 


Re s 

U c 

d 

V 

t X X{ y X 4 

N x xN y xN z 

N P 

100 

1 

0.8 

8 • 10~ 3 

2nx2nx I6n 

256 x 256 x 2048 

« 10 6 

400 

1 

0.8 

2 • 10- 3 

2nx2nx 16 n 

256 x 256 x 2048 

* 10 6 

1000 

1 

0.65 

6.5 • 10- 4 

2nx2nx 16 n 

512x512x4096 

* 10 6 


TABLE I. Parameters of the numerical simulations. Re s = U c d/v : Reynolds number of the obstacle, U c : inflow velocity, d: 
diameter of the spherical obstacle, v: kinematic viscosity, l x , l y and € z : edge lengths of the computational domain (z being in 
the direction of the mean flow), N x , N y , and N z : number of collocation points in the respective directions, N p : number of small 
particles. 

The simulations are set up in the following way. First, we integrate the Navier-Stokes equations 0 with a sphere 
located at Xs = (tt, tt, 2 n) while imposing a homogeneous inflow with U c = 1 until a (statistically) stationary state is 
reached. We then inject a stream of inertial particles by placing randomly many particles per time step at the entrance 
of the domain. They are initialized with the velocity Uc = 1 of the fluid. During the simulation we remove all the 
particles that are touching the sphere or reaching the end of the computational domain. On average the domain is 
filled with approximately one million particles. The fluid and particle trajectories are then stored and analyzed during 
four times the time needed for the fluid to travel across the numerical domain. 


III. DENSITY INHOMOGENEITIES IN THE WAKE 
A. Average concentrations and voids 

Small inertial particles that pass over the spherical obstacle are influenced by both its boundary layer and its wake. 
The rotating structures in the wake eject particles and create voids and concentrations (see again Fig. [I]). The mean 
(temporal average) particle number density is shown in Fig. [2] for the three different Reynolds numbers we have 
considered. Due to the axial symmetry of the problem, we adopt a cylindrical coordinate system (r, z, 0) (where z is in 
the direction of the inflow Uc) and average over the angle 0. One observes regions where the particle concentration 
is either significantly over or under its value in the inflow stream. 

Let us outline the main findings: A cylindrical region where small particles concentrate is emerging at the sur¬ 
face of the obstacle while a shadow with a very low concentration extends downstream. In the case of the two 
highest Reynolds numbers (turbulent wakes), this shadow region is drastically reduced and, moreover, another over¬ 
concentration is created further downstream. This second high-concentration region around the axis of symmetry 
is accompanied by a conic region of reduced concentration spreading outwards. Both the high as well as the low 
concentrations are more pronounced for the higher Reynolds number. For St = 3.9 and Re s = 1000, we observe for 
example an over-concentration in the wake that is approximately 1.5 times larger than the mean concentration. 


B. Stream-wise and radial profiles 

The spatial dependence of these concentrations can be seen with more details in Fig. [3] (left). Particles with a 
small Stokes number St only create a short shadow (a few times the obstacle diameter d for St = 0.12) where the 
concentration is significantly depleted. Rapidly, at only several diameters further downstream, the density exceeds the 
inflow density and reaches a maximum. Beyond this point, the over-concentration slowly reduces. This is because t T] 
increases (and therefore St^ = r/r^ reduces) as a function of the distance from the sphere so that small-5f particles 







4 



FIG. 2. Mean particle concentration for St = 0.4 and Re s = 100 (top), St = 0.4 and Re s = 400 (middle) and St = 0.49 and 
Re s = 1000 (bottom). 
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FIG. 3. Left: Mean particle concentration on the axis of symmetry as a function of the distance from the spherical obstacle for 
several St and Re s = 1000. Data has been smoothed by a sliding window average with a window size approximately equal to d. 
Right: Position in the wake where the small particle density recovers the inflow density, corresponding to the transition between 
the shadowed region and the over-concentration. 


behave more and more like tracers and stick to fluid elements. They are transported and mixed by the wake turbulent 
fluctuations and this evens out possible concentrations. 

The density profile for particles with a larger Stokes number St is similar but shifted further downstream due to 
higher inertia. The concentrations might persist in the far wake as turbulence decays downstream so that heavy 
particles do not diffuse anymore and their density might just be advected by a quasi-uniform flow. Also, it seems 
that the maximum of concentration increases with St. However, these two observations cannot be unquestionably 
confirmed from our simulations, given the limited span in the stream wise direction. A more precise investigation of 
these aspects would require a much longer computational domain and is kept for a future study. 

The starting point of the over-concentration region (where the particle density exceed that of the inflow) is shown in 
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Fig. [3] (right) as a function of St. It follows a power-law with an exponent which is close to that measured by Jacober 
and Matteson. 11 The decrease of the over-concentration as a function of the distance from the obstacle is reported in 
Fig. [4] Here, only small-SY particle are considered as the position of the maximal density lies outside the domain for 
large St. The decay follows a power-law with an exponent close to -2/3. 



FIG. 4. Left: Mean particle concentration deficit on the axis of symmetry as a function of the distance from the obstacle center 
for several St. Data has been smoothed by a sliding window average with a window size of approximately d. Data is shown for 
Re s = 1000. The black line corresponds to a decay oc z~ 2/3 . 




FIG. 5. Radial density profiles for various distances downstream the spherical obstacle ( Re s = 1000) and for St = 0.12 (left) and 
St = 3.94 (right) 

The radial profiles of the particle density are shown in Fig. [5] for two different Stokes numbers. At small distances 
z from the obstacle, the high-concentration jet emerges at a radial distance r « d from the axis of symmetry. Further 
downstream, the transition from regions of reduced to increased particle concentrations are also situated at r « d. 


C. Concentration mechanisms 

For all Reynolds numbers of the spherical obstacle, a high-concentration cylindrical jet of particles emerges directly 
from the surface of the sphere. The reason for this is the converging flow in front of the obstacle (see the sketch in 
Fig.®. Particle trajectories that start upstream at locations close to the axis of symmetry (small r’s) come close to 
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those starting with a larger r. This is due to reduced accelerations of the inertial particles compared to fluid elements. 
The observed number densities in the jet region decrease with increasing Stokes number St. This is because the 
collision rate of small particles with the sphere increases with St. The particles which have collided are removed from 
the flow and are thus missing downstream and in the jet region. 



FIG. 6. Sketch showing the mechanisms leading to over- and under-concentrations of small particles in the wake of a spherical 
obstacle (red half-disk). The upper half shows the situation for laminar wakes while the lower part relates to turbulent wakes. 
Streamlines are shown in black. In the case of a turbulent wake, the stream lines of the mean flow are shown as dashed lines. Blue 
lines correspond to trajectories of heavy particles. In the turbulent wake particle positions are indicated by blue points. Green 
spirals sketch the position and rotation of detached vortices. 

The low-concentration shadow behind the obstacle originates from the fact that the dynamics of inertial particles 
lies behind that of the fluid. Fluid streamlines converge downstream the recirculation region but inertial particles 
overshoot. The resulting shadow persists far away in the wake (in the laminar case), because of the turbulent decay 
downstream and the associated decrease of the fluid radial acceleration. Particles get less and less driven toward the 
center axis of symmetry. 

How can one understand the emergence of the over-concentration regions in the far turbulent wakes? Jacober and 
Matteson 11 suggested that they originate from the jet region that is drawn around the sphere and collapses behind it. 
However, if this was the case, a laminar wake would also produce similar concentrations and we do not observe this 
effect. The mean flow in the turbulent wake is not responsible for these concentrations either. We indeed integrated 
trajectories of inertial particles according to but transported by the mean velocity field u(x), i.e. the temporally 
averaged velocity field of the Re s = 400 simulation (not shown). Up to several diameters downstream from the fixed 
sphere, the concentration profile is very similar to the time dependent simulations. Notably, the jet is reproduced. 
However, this approach does not produce the high-concentration region in the far wake. 



FIG. 7. Mean particle concentration in pseudo colors for St = 1.6 (r = 1.28) particles and Re s = 400. Position of the maximal A 
as a function of the downstream distance z is shown as a dashed line. Solid black lines indicate contours of the local Kolmogorov 
time-scale r^. 

We suggest that the concentration mechanism is provided by the rotating structures in the turbulent wake. Yao et 
tf/Pi already observed instantaneous concentrations in a two dimensional flow past a disk and found that particles con¬ 
centrate between periodically shaded vortices. Tang et alW and Yang et alW^ studied numerically and experimentally 
the influence of advected large scale structure on particle dispersion in plane mixing layers and wakes. They found 
that the development of highly organized patterns of inertial particles is determined by the evolution and interaction 
of these structures. Via fractal correlation dimension measurements they showed that particles with St = 1 are the 
most influenced. 

In order to study with more details the effects of the wake rotating structures on particle concentrations, we com- 
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puted the temporally averaged Okubo-Weiss parameter A (defined in the caption of Fig. [TJ. It is a measure of the 
strength of rotation vs. stretching in the fluid velocity field. A is positive in regions where the velocity gradients are 
dominated by vorticity and negative in those where strain is dominant. In a turbulent flow, heavy inertial particles 
tend to move from regions where A is positive to those where it is negative. 7 The maximum of A as a function of the 
streamwise position is shown as the dashed line in Fig. [7] One clearly observes that it correlates well with the region 
of under-concentration. 

Another observation is that, when the particle inertia is increased, the starting point of the concentration region is 
shifted downstream (compare Fig. [3]). This change can be explained by the fact that particle ejection from rotating 
regions is the most effective for particles with a Stokes number St r] = r/r^ around unity. The local Kolmogorov 
time scale increases downstream (see solid black contour lines in Fig. [7]). Hence, the position where St ^ l shifts 
downstream when particle inertia increases. 


IV. INTER-PARTICLE COLLISIONS 

The goal of this section is to qualitatively estimate the occurrence and position of possible collisions between 
particles with equal inertia. We will limit the discussion to the flow Re s = 400 and we focus on mean quantities, 
leaving the influence of fluctuations for a future work. Based on the mean pair density and on the mean approaching 
velocity between particles with the same Stokes number, one can estimate the collision rate is proportional to the 
product 


Pc(z, r) = -<n_(z, r ) 2 > < W-(z , r)), (4) 

where n 2 is the coarse-grained pair density of approaching particles and w_ the associated mean longitudinal velocity 
difference. For both n 2 and w_, only particles with negative longitudinal velocity differences are taken into account. 
Such an estimate for collision rates relies on the ghost-particle approximation^ where possible particle encounters 
are counted but not performed. It relies on providing an arbitrary radius to the particles that we have here chosen to 
be 0.2 d. All the coarse-grained quantities used in 0 are thus defined for that scale. We observe from our simulations 
that the pair density ( n 2 _) has the same structure as the mean single particle concentrations (not shown): A high 
concentration jet region in the vicinity of the large sphere for particles with a small inertia and a high concentration 
region downstream, more pronounced for particles with a large inertia. 



z/d 



FIG. 8. Top: Average velocity difference (w_(z, r)) of approaching particles with St = 3.2 , coarse-grained at a scale 0.2 d. Bottom: 
Estimate p c , defined in & for binary collisions of particles with the same St = 3.2. It is represented in logarithmic units in order 
to emphasize its variation rather than only its maximum in the jet region. 

The average longitudinal velocity differences are shown in Fig. [8] (top). The highest approaching velocities are 
located in the jet region, where particle trajectories are compressed. In the wake region where vortices detach, one 
still observes some non-vanishing values of the approaching velocities, which are much less strong than in the jet 
region. Further in the wake, the velocity differences become weak due to the decay of turbulent fluctuations. 

The estimate p c for collision rates is shown in Fig. [8] (bottom). The jet region provides the highest probability 
for particle encounters. There, both the particle density as well as the particle approaching speed are enhanced. 
With increasing distance from the obstacle the particle collision region moves towards the axis of symmetry and the 
probability of collisions decreases. The dependence of the collision probability as a function of the distance from 
the obstacle can be seen more clearly in Fig. [9] where we show the sum of p c over planes (r, 6) as a function of the 
streamwise direction z. One clearly observes that the particles with a small inertia collide preferentially in the jet 
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region, in the vicinity of the sphere. The collisions between particles with a larger inertia are rather enhanced in the 
wake. We indeed see in Fig. [9] that non-vanishing values of p c persist at downstream distances as far as tens of the 
obstacle diameter. 



FIG. 9. Sum of p c over planes (r, 6) in directions perpendicular to the stream, as a function of the downstream distance to the 
obstacle. Three different strengths of particle inertia are shown, as labeled. 


V. SUMMARY AND OUTLOOK 

In this work we have analyzed how a stream of inertial particles is influenced by the boundary layer and the wake of 
a spherical obstacle. We make use of the results of three direct numerical simulations of the incompressible Navier- 
Stokes equations using a combination of a Fourier-spectral method and a penalty technique. We considered three 
different Reynolds numbers for the sphere, Re s = 100, Re s = 400, and Re s = 1000. The average concentration of 
small inertial particles shows different characteristics in the wake of the spherical obstacle, depending on the particle 
Stokes number St and on the sphere Reynolds number Re s : For all Re s we find that particles with weak or moderate 
inertia cluster at the outer edge of the obstacle creating a narrow cylindrical jet of particles. Behind the obstacle, inertia 
is responsible for the formation of a shadow where particle concentration is strongly reduced. When the wake of the 
obstacle is laminar, these over- and under-concentrations persist far downstream while they are rapidly dissipated by 
fluctuations for a turbulent wake. Moreover, we find that in the turbulent case, there is a region appearing a couple of 
diameters downstream the sphere where the particle concentration exceeds that in the inflow. This can be explained 
by resonant ejections of particles from detached and advected vortices with a local turn-over time of the order of 
the particle response time r. 

Estimates of the probability of collisions between identical particles show that particles with a weak inertia tend to 
collide in the jet region near the sphere while those with a strong inertia collide preferentially further downstream. In 
a future work it will be interesting to examine such collision statistics in the poly disperse case, i.e. between particles 
with different sizes and response times. We expect that in the turbulent wake the relative velocity of particles with 
different inertia will be larger than for particles of the same type. We also plan to take coagulation or coalescences 
into account and analyze how the particle size distribution is affected by the wake of the spherical obstacle and how it 
depends on the Reynolds number of the latter and on the downstream distance. This can have important applications, 
for instance in planet formation. It seems well founded to assume that large bodies, such as planetesimals or planetary 
embryos, will concentrate dust in their wake and will thus foster accretion and the formation of other objects. 
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